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In recent works, the idea of time-reversed laser oscillation has been proposed and demonstrated 
to realize a two-channel coherent perfect absorber [Y.D. Chong et al., Phys. Rev. Lett. 105, 
053901 (2010); W.Wan et al., Science 331, 889 (2011)]. Here the time reversal of optical parametric 
oscillation in a nonlinear x'^' medium is considered and shown to realize a coherent perfect absorber 
for colored incident signal and idler fields. A detailed analysis is presented for the time-reversed 
process of mirrorless optical parametric oscillation in the full nonlinear regime. 
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Introduction. Time-reversal symmetry plays an impor- 
tant role in classical electromagnetic theory and nonrela- 
tivistic quantum mechanics [T] . Such a symmetry implies 
that, if a particular physical process is allowed, then there 
also exists a time-reversed process, which is obtained by 
changing the sign of the time variable in the dynami- 
cal equations. In recent works, the time-reversed process 
of laser oscillation has been explored, both theoretically 
[2^ and experimentally |3], and shown to realize a coher- 
ent perfect absorber (CPA). Under additional symmetry 
constraints, the processes of lasing and CPA can even 
occur simultaneously, and a CPA-laser device can be in 
principle realized [4i |5j . Since in a steady-state process 
time reversal corresponds to interchanging incoming and 
outgoing fields, the time-reversed process of a laser insta- 
bility corresponds to perfect absorption of incoming light 
fields. Such a phenomenon is rather general and expected 
to occur for the time reversal of other optical instabili- 
ties. Among these, optical parametric oscillation (OPO) 
provides an important example of instability in nonlin- 
ear optics, which is at the core of the operation of such 
important devices as optical parametric oscillators [BHS] ■ 
In a parametric oscillator at the onset of instability, a 
pump field at frequency uj^ propagating in a nonlinear 
T^(2) crystal generates down-converted fields at frequen- 
cies wi (the signal field) and uj2 = uj^—oji (the idler field), 
provided that the phase matching (momentum conserva- 
tion) condition is satisfied. Generally, at least one of the 
signal or idler waves resonates in an optical cavity. In the 
time-reversed process, exactly the opposite occurs: two 
coherent signal and idler fields, with appropriate ampli- 
tude and phase relationship incident onto the pumped 
nonlinear crystal, are completely annihilated, thus re- 
alizing a kind of "colored" CPA device. While a CPA 
based on time reversal of a laser usually converts the 
absorbed light into heat [H [3], in the time reversal of 
OPO the annihilated photons are converted into light of 
different frequency (the sum-frequency field). Moreover, 
since the OPO instability and its time-reversal process 
are allowed in the same medium, an ordinary OPO de- 
vice can behave at the instability like a CPA-laser device 
01 E] , i.e. it can simultaneously emit outgoing colored 
fields and absorb incoming colored fields with appropriate 
phase relationships. This important property shows that 



a conventional OPO system at threshold is a CPA-laser 
of a certain kind. In this Letter we consider specifically 
the time-reversed process of OPO for counterpropagat- 
ing photons, i.e. the time-reversed process of backward 
OPO [H Uni [11] , which enables a simple analytical study. 
Backward OPO describes an interesting three-wave in- 
teraction process in a quadratic medium which has the 
unique property of automatically establishing distributed 
feedback, thus allowing mirrorless OPO. Mirrorless OPO 
was theoretically predicted in the early days of nonlin- 
ear optics pp] , however its experimental demonstration 
has been reported only in recent years using short-period 
quasi-phase-matching (QPM) nonlinear crystals 12, 13]. 
Backward parametric interaction is also useful for the 
generation of entangled twin photons [13] . 

Time-reversal of OPO instability: general aspects. We 
consider interaction of three quasi-monochromatic waves 
at carrier frequencies wi (signal field), ui2 (idler field) 
and a;3 — ioi + (jJ2 (pump field) in a nonlinear x^"^^ 
medium in the framework of classical nonlinear optics 
equations. In the scalar approximation, the electric field 
£[r,t) satisfies the wave equation V^f — {l/cQ)d^£ — 
{1 / eQC^)dfV , where cq is the speed of light in vac- 
uum, eg is the vacuum permittivity, is the spa- 
tial Laplacian, and V — Vl + Tnl is the polariza- 
tion of the medium, which includes the linear {Vl) and 
nonlinear (Vnl) terms. For nonlinear interaction in a 
quadratic medium with QPM, one can assume Vnl ~ 
eox^^^S^, where x^^H^") is the spatially-modulated 
second-order susceptibility of the medium. After setting 
£{r,t) = {l/2)[Ei{r,t)exp{iujit) + E2{r,t) exp{iuj2t) 
E3{r, t) exp(icj3t) -l-cc], where £';(r, t) varies slowly with 
respect to time t over one oscillation cycle l-Kjioi (I = 
1, 2, 3), neglecting group velocity dispersion and indicat- 
ing by e; (r) the relative dielectric constant of the medium 
at the frequency , the envelopes Ei satisfy the coupled 
wave equations (see, for instance, |15l ) 

[ V? + ^ - ^^.'l E,.2 ^ ~^\x'^El,E, 

^l^^^-^^d^E, = -%^-^E,E2 (1) 
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where ki — uini/ca, ni and vi {I = 1,2,3) are the wave 
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FIG. 1. (Color online) Schematic of (a) an OPO instability, 
and (b) its time-reversed process, corresponding to a colored 
CPA. 



number, bulk refractive index and group velocity at fre- 
quency uji, respectively. For a lossless medium, the di- 
electric constant e; (r) is real- valued and can be generally 



written as e;(r) 



2niAn;(r), where Ani(r) is a 



small correction to the bulk refractive index n/ that ac- 
counts for e.g. index guiding or cavity effects for the sig- 
nal and/or idler fields, like in singly- or doubly-resonant 
OPOs. Equations (1) should be supplemented with the 
boundary conditions at the crystal (or cavity) bound- 
aries. In an OPO geometry, a traveling-wave pump field 
shines on the crystal, and the signal and idler fields are 
spontaneously generated at the OPO instability as sta- 
tionary outgoing waves. Time reversal of the underlying 
equations implies that, if Ei(r) is a stationary solution to 
Eqs.(l) at the OPO instability point, then ^^^(r) is still 
a solution to the same equations. Since the complex con- 
jugation replaces outgoing waves with incoming waves, 
the time-reversed process of OPO is the perfect annihila- 
tion of the incoming signal and idler fields, i.e. a colored 
CPA is realized. A schematic of the OPO instability 
and of its time-reversal is depicted in Fig.l. It should 
be noted that the present analysis is strictly valid only 
for classical fields, and phenomena involving few photons 
as well as quantum noise are not considered in such a 
semiclassical approach. If quantum fields are considered, 
spontaneous parametric down-conversion occurs, and the 
complete annihilation (conversion to the pump field) of 
both incident fields cannot be strictly reached. 

Time-reversed backward OPO. In the following, we will 
focus our analysis to the backward OPO configuration in 
a QPM crystal of length L, which allows for mirror less 
OPO [i.e. An/ = in Eqs.(l)]. In the plane-wave ap- 
proximation [TT] , for a backward configuration the pump 
and signal fields copropagate along the positive z direc- 
tion, whereas the idler wave is counter-propagating along 
the negative z direction, as shown in Fig. 2 (a). Phase 
matching of the interaction is ensured by a QPM grating 
X^^)(z + A) = x'^^K^) with period A satisfying the first- 
order QPM condition 2tt/A= |A;3-|-fc2 — ^il- After setting 
Ei^3{z,t) = (l/K0(wi,3/"-i,3)^''^ai,3(a;,t)exp(-ifci^3z), 
E2{z,t) — {l/Kl){u}2/n2)^^'^a2{z,t)exp{ik2z) for the 
forward- and backward-propagating waves and introduc- 
ing the normalized spatial propagation distance Z — z/L, 
from Eqs.(l) the following coupled-mode equations can 
be obtained for the dimensionless envelopes ai after av- 
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FIG. 2. (Color online) (a) Schematic of backward parametric 
interaction of pump, signal and idler waves in a nonlinear 
QPM crystal, (b) Behavior of the parametric gain $ versus 
the phase difference of incident signal and idler fields for 
increasing values of the normalized input pump amplitude 
aao in the no-pump-depletion approximation and below the 
OPO threshold. Curve (1): ago ~ 7r/6; curve (2): aso — 7r/3; 
curve (3): aao — n/2.1. The dip at tf) — n/2 is the signature 
of colored CPA. 



eraging the rapidly-oscillating terms in z |TT1 [TS] 



± dzai^2 + {L/vi,2)dtai^2 = -ia^^a^i 



(2) 



where k = (l/2co)(wia;2W3/riin2"-3)^^^XAf and xm is 
the Fourier amplitude of x {^) at the spatial frequency 
27r/A. Without loss of generality, k is assumed to be real- 
valued and positive. The intensities /; of the three waves 
are given by/; = eoConin2n3XiX2X3\ai\'^ /[Zn'^ L'^xl-iniXi], 
where A; — 2t:cq / loi are the (vacuum) wavelengths of the 
three waves [l — 1,2,3). To investigate the time rever- 
sal of OPO, we assume that a continuous-wave pump 
with amplitude 03(0) = 039 is injected into the nonlin- 
ear crystal, together with two signal and idler waves with 
amplitudes ai(0) and 02(1), respectively. 
Let us first assume that the amplitudes of injected signal 
and idler fields are much smaller than the pump ampli- 
tude 030 (small-signal analysis). Then in the no-pump- 
depletion approximation, a3(Z) ~ aso, Eqs.(2) reduce 
to a system of two coupled linear equations for the am- 
plitudes ai{Z) and a^iZ) jH|, which yield the algebraic 
relation (ai(l), 03(1))^ = 7W(ai(0), 03(0))^ between the 
amplitudes of signal and idler fields at input {Z = 0) and 
output {Z = I) crystal planes. Assuming without loss 
of generality a positive and real- valued amplitude 039, 
the coefficients of the transfer matrix Ai read explicitly 
Mil = M22 = cos(a3o) and A^i2 = A^2i = -«sin(a3o). 
The vanishing of A^22 as the pump amplitude is increased 
from zero defines the mirrorless OPO threshold, which is 
attained at 030 — 7r/2 (see, e.g., [51|TT]). It is worth ob- 
serving that the transfer matrix M has the same struc- 
ture as that encountered in the PT-symmetric CPA-laser 
device, recently proposed in Refs.[31 [S]. In particular, 
since Mn = M22, the simultaneous vanishing of Mu 
and A^22 at the OPO instability threshold indicates that 
the system can simultaneously emits coherent signal and 
idler waves and acts as a colored CPA for appropriate 
excitation. This means that an OPO at the instability 
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threshold basically realizes a kind of colored CPA-laser 
system. Like in the CPA-laser device, colored CPA re- 
quires coherent excitation of the crystal with signal and 
idler fields of appropriate amplitude and phase relation- 
ship, and the degree of absorption/amplification can be 
conveniently tuned by varying the relative phase of the 
two incoming beams. In fact, let us consider a pump 
level close to (but below) the threshold for OPO insta- 
bility, and let us assume that the crystal is simultane- 
ously illuminated, from opposite sides, by signal and idler 
waves with the same amplitude A and phase difference (/>, 
namely let us assume ai(0) = Aexp^icj)) and 02(1) = ^■ 
The parametric gain for signal and idler waves after prop- 
agation in the crystal is given by — \ai{l)/A\'^ ~ 
\a*{Q)/A\^ = {cos2(0) + [sin(0) - sin(a3o)]n/ cos2(a3o). 
A typical behavior of $((/') versus the phase difference 
is shown in Fig. 2(b) for increasing values of the pump 
level a^Q. Note that, at (j> — tt/2 a, dip is observed in 
the gain curve, with $(7r/2) — > as the threshold of 
OPO is approached. Such a dip is the clear signature 
of CPA. Far from (j) — 7r/2, the system behaves con- 
versely like a parametric amplifier ($ > 1), with a para- 
metric gain that diverges as the OPO threshold is ap- 
proached. The curves shown in Fig. 2(b) clearly indicate 
that an optical parametric amplifier is capable of attenu- 
ating strongly as well when coherently illuminated from 
both sides, for an appropriate relative phase of the in- 
coming fields. Such a result basically stems from the 
well-known phase-sensitivity of the frequency conversion 
process, which controls the flow of energy from the pump 
to the signal and idler waves. 

The previous considerations are based on a small-signal 
analysis of Eqs.(2), which neglects pump depletion ef- 
fects. However, a more comprehensive analysis should 
consider the full nonlinear regime, i.e. should account 
for pump depletion effects. Such an analysis is of ma- 
jor importance to check the stability of the CPA regime 
against the spontaneous emission of outgoing signal and 
idler waves. The exact solution corresponding to CPA in 
the nonlinear regime can be constructed following a sim- 
ilar procedure used in the theory of backward OPO [TTj, 
but with modified boundary conditions. Let us search 
for a stationary solution to Eqs.(2) of the form 

di{Z) = aoSin[6'(Z)] , a2{Z) = -mo cos[6i(Z)], (3) 

where qq is a complex parameter and the angle 0{Z) is 
an unknown function. Substitution of Eq.(3) into Eqs.(2) 
yields the following coupled equations for a-i{Z) and 9{Z) 

dz0 = a3, 9za3 = -(|ao|V2)sin(20). (4) 

The colored CPA regime is attained by imposing the 
boundary conditions ai(I) — 0,2(0) — for the signal 
and idler fields, and 03(0) = 030 for the pump field. The 
former conditions can be satisfied by taking 0(0) — 7r/2 
and 9(1) = mr, where n = 1,2,3,.... Note that in this 
case |ao| corresponds to the normalized amplitudes of the 
counter-propagating signal and idler fields incident onto 



the crystal at the Z = and Z = 1 facets, respectively. 
For an assigned pump level 030, the value of |ao| has 
to be consistently determined from Eqs.(4) with the ap- 
propriate boundary conditions. In fact, from Eqs.(4) it 
follows that the angle 0{Z) satisfies the pendulum equa- 
tion {(Pe/dZ'^) + (|ao|V2)sin(26') = 0, which admits of 
the constant of motion (energy) E = (l/2){dd/dZ)'^ — 
(|aop/4) cos(20). Hence the following relation holds 

/ , = / dZ = l. (5) 

Je{o)=^/2 \/2£;-h (|aol72) cos(26i) Jo 

The value of the energy E is determined by the boundary 
condition at ^ = and reads E = (I/2)a|o + |aoP/4. Af- 
ter setting 77 = |ao/'^30p5 from Eq.(5) one readily obtains 



(v/TT^) a3o = (2n - l)K YTv) 

where K{£,) = J^^^ dq/[l — sm^{q)]^/'^ is the complete 
elliptic integral of the first kind. For a given order n, 
Eq.(6) implicitly defines the relation between |ao| and 
030, which is depicted in Fig. 3(a) for the lowest-order 
(n = I) state. Hence, to realize the CPA condition, two 
conditions should be satisfied: (i)the relative phase be- 
tween signal and conjugate idler waves is fixed accord- 
ing to Eq.(3); (ii) for a given injected pump level 030, 
the normalized signal and idler wave amplitudes injected 
into the crystal from the two sides should take the same 
value |ao| defined by the curve of Fig. 3(a). In partic- 
ular, for small values of |ao| the amplitude 030 of the 
pump wave that realizes the CPA condition is exactly 
the threshold value 7r/2 of OPO instability, according to 
the small-signal analysis previously discussed. The ex- 
plicit behaviors of ai{Z) along the crystal are then ob- 
tained by solving Eqs.(4) and using Eqs.(3). Note that 
the phase of the amplitude Oo , which defines the absolute 
phases of the injected signal and idler fields , remains 
fully undetermined, so that if {ai{Z),a2{Z),a3{Z)} is 
a stationary solution to Eqs.(2) corresponding to CPA, 
then {ai(Z) exp(iip), a2{Z) exp(— i(^), a^{Z)} is a solution 
as well for an arbitrary phase ip. As an example. Fig. 3(b) 
shows the numerically-computed behavior of the normal- 
ized field intensities |a;(Z)p for pump, signal and idler 
fields along the crystal for the lowest-order CPA solution 
(n = 1) corresponding to = 5 [point C of Fig. 3(a)]. 
Note that the annihilation of the signal and idler photons 
corresponds to the creation of pump photons, rather than 
to the depletion as in the backward OPO [II . Higher- 
order {n > 2) CPA solutions show oscillating power ex- 
change among signal, idler and pump fields along the 
crystal. As an example. Fig. 3(c) shows the behavior of 
|a;(Z)p along the crystal for the n = 2 mode correspond- 
ing to ?7 = 5. 

The stability of the CPA solution can be investigated 
by a standard linear stability analysis [16]. After setting 
oi{Z, t) = a{Z) + Soi{Z, t) in Eqs.(2), where 5ai{Z, t) are 
small perturbations [l = 1,2,3), the most general solu- 
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FIG. 3. (Color online) (a) Behavior of the input pump am- 
plitude aso versus the amplitude |ao| of input signal and idler 
fields corresponding to the lowest-order CPA solution, numer- 
ically computed from Eq.(6) with n = 1. The points A, B and 
C in the figure correspond tor; = 0.5,r7 = l and = 5, respec- 
tively, where r\ = |ao/a3o|^- (b) Behavior of the normalized 
intensities |a;(Z)|^ of pump, signal and idler fields inside the 
nonlinear crystal for the lowest-order [n — 1) CPA solution 
corresponding to 77 = 5 [point C of Fig. 3(a)]. (c) Same as (b), 
but for the second-order (n = 2) CPA solution with = 5. 
(d-f) Loci of the eigenvalues A of the linearized equations for 
the lowest-order CPA solution corresponding to points A, B 
and C of Fig. 3(a). (g) Same as (d-f), but for the second-order 
CPA solution of Fig. 3(c). In the computation of the eigen- 
values, the same group velocity vi — V2 = V3 = v has been 
assumed for the three waves. 



tion to the linearized equations 5ai {Z, t) is given by an ar- 



bitrary superposition of solutions of the form 6ai{Z, t) = 
ui{Z) exp(Ai) -I- vl{Z) exp(A*i), where the vector of per- 
turbation amplitudes v = (ui, U2, U3, Vi,V2, v^Y' satisfies 
the linear system {dv/dZ) = t/(Z, A)v and A is a com- 
plex parameter whose real part determines the growth or 
decay rate of the perturbation. The explicit form of the 
6x6 matrix Q{Z, A) is given in the supplementary ma- 
terial |16) . The eigenvalues A are then determined after 
imposing the boundary conditions for the perturbations 
at the Z = Q and Z = 1 planes [16]. Owing to the 
phase invariance of the CPA solution, a neutrally stable 
mode with eigenvalue A = is always found. An insta- 
bility arises provided that the real part of at least one 
of the other eigenvalues is positive. Extended numerical 
simulations indicate that the lowest-order (n = 1) CPA 
solution is stable at any pump level, whereas higher-order 
(n > 2) CPA solutions are linearly unstable. As an ex- 
ample, Figs.3(d-f) show the numerically-computed loci 
of the eigenvalues A for the lowest-order (n — 1) CPA 
solution at ?7 = 0.5, rj ~ 1 and 77 = 5, respectively, corre- 
sponding to points A, B and C in Fig. 3(a). Note that all 
the eigenvalues lie in the lower-half complex plane, and 
thus there are not unstable modes. Conversely, Fig. 3(g) 
depicts the eigenvalues for the n = 2 CPA solution at 
77 = 5, indicating the existence of an unstable mode as- 
sociated to the eignevalue with a positive real part. 

Conclusion. In this work, time reversal of optical para- 
metric oscillation has been proposed and shown to realize 
a kind of "colored" coherent perfect absorber. A detailed 
analysis has been presented for the time reversal of mir- 
rorless OPO. At threshold and in the no-pump-depletion 
limit, it has been shown that the OPO can simulta- 
neously emit outgoing coherent signal and idler waves, 
and perfectly absorb incoming coherent signal and idler 
waves, thus realizing a kind of CPA-laser device [HIS]. As 
compared to previous studies on CPA and CPA-laser sys- 
tems [J'-^, CPA operation has been extended here into 
the nonlinear regime. In particular, for the mirrorless 
OPO a proof of the stability of the colored CPA opera- 
tion for the fundamental CPA mode has been presented. 
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